OPERATIONAL METHOD FOR FINITE DIFFERENCE EQUATIONS 
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Abstract. In this article I present a fast and direct method for solving several types of linear 
finite difference equations (FDE) with constant coefficients. The method is based on a polynomial 
form of the translation operator and its inverse, and can be used to find the particular solution of the 
FDE. This work raises the possibility of developing new ways to expand the scope of the operational 
methods. 
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1. Finite Difference Equations (FDEs). These equations describe the re- 
lationship between the present value of a function and a discrete set of n previous 
values: 

f{t)J{t + l),...J{t + n) 

The solution of an FDE is a function /(i), with t E Z, that satisfies the equation. 

Typically, the n known values for a FDE of degree n are referred to as initial 
conditions. In order to establish a unique solution, it is necessary to know the values 
of the initial conditions as precisely as possible. They are habitually defined at equal 
intervals starting from t = 0: /(O), /(I), /(2), . . . , /(n). In this notation, the FDE 
predicts f{n). 

1.1. Linear FDEs. A linear FDE can be expressed as follows: 

ao{t)y{t + n) + ai{t)y{t + n - I) + ■ ■ ■ + an-i{t)y{t + 1) + a„(t)?/(t) = 

To solve a linear equation with constant coefficients (LECC), 

aoyit + n) + aiy(t + n - 1) H h a„_iy(< + 1) + a„?;(t) = (/)(t) , 

we define the Translation Operator T: 

y{t + 1) - Ty{t); y{t + k) ^ T{T{- ■ ■ {T y{t)) ■■■))= T^y{t) ^ 

fe 

aoT"2/(t) + aiT"-iy(t) + • ■ • + an-iTy{t) + a„y(t) = ^ P{T)y{t) = 
and look for a general solution of the form 
(1.1) 

yGeneral\J ) — V Homogenous [J ) i" V Particular {t)^yG{t)^yH{t)+yp{t) . 

Several methods for solving LECCs exist. In this paper, I develop an operational 
method to find a particular solution. 

2. Operational Method for Finite Difference Equations. The basic idea 
of the operational method is simple. If we wish to define the polynomial operator 
P{T), is is possible to establish its properties? We begin with the following axiom: 



(2.1) P(T)g(t) = f (t) ^ (t) = g(t) 
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2.1. Properties of 
1. Linearity 



P(T) 



P(T) am+pgit) 



ao af{t + n)+Pg{t + n) +ai af{t + n-l) + l3g{t + n-l) + 



+ ■■■ + an-i af{t + 1) + pg{t + 1) +a„ af{t) + pg{t) = 



aof{t + n) + aif{t + n - 1) + • • • + a„_i/(t + 1) + a„/(i) 



+/3 ao5(t + n) + aiS((f + n - 1) H 1- an-ig{t + 1) + a„g(t) = 



= aP(r)/(t)+/3P(r)5(0 



Therefore: 



P(T) af{t) + pg{t) = aP{T)f{t) + pP{T)g{t) ^ 



(2.2) 



P(T 

2. Inverse Translation 



rf(t)+^g(t) 



--p^f(t) + /3p^g(t) 



T"/(t) = /(t + n)^/(i) = ^/(i + n) 



We change t ^ t — n: 
(2.3) 



-f(t) = f(t-n) 



3. Unity 

Let y{t) = t, with n £ N. The expression y(f + 1) — y(t) is equivalent to 
f+l-i = 1, so the FDE y{t + l)-y{t) = 1 has the solution y{t) = t. Finally, 



(2.4) 



T- 1 



(l)=t 



4. Propagation 

Let us solve the two simplest cases directly: 



y{t + 1) - y{t) = t 



t _tt-l 
T-1 ~ 1 ' 2 
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y{t + 2) - 2y{t + 1) + y{t) = t j^^^ = I " ^ ' ^ " 
Generalizing the progression, we obtain 

t t t-l t-2 t-n t-i 



Therefore, 



(T-1)" 1 2 3 n + 1 -•-•■« + 



t - i 



^ i=0 

2.2. Equations where (?i(t) = A*. Consider the polynomial 
P(T) = aoT" + aiT"-i + • • • + a„_iT + a„ . 

We have 

P(T)A* = aoT"A* + oiT"-1A* + • ■ • + a„_iTA* + o„A* = 
= aoA*+" + aiA*+"-i + • • • + a„_iA*+i + a„A* = 
= A*(aoA" + aiA"-i + ■ • ■ + a„_iA + a„) = A*P(A) 

Therefore, given that -P(A) 7^ always holds, we consider the property 

(2.6) p(T)A* = A*P(A) 4^ -1-A* - 



P(T) P(A) 

Example 2.1. Find a particular solution yp{t) of the finite difference equation 
y{t + 2)-5y{t + l)+4y{t) = 3* . 

Solution: 

y{t + 2) - 5y(i + 1) + 4y{t) = 3* ^ (T^ - 5T + 4)y(i) = 3* 

Replacing: 

= r^-5r + / = 3^-5-3 + 4 ^ = -| 

□ 

2.3. Equations where (j){t) = cosuTrt or (/)(t) = sinuTrt. We begin with Euler's 
formula: 

gWTTtl _ _|_ ^ gin j^TT^ _ 

In particular, we consider n e N: 

e"'^' = cosnTT + isinnTT = (-1)" + i • = (-1)" 

We apply the operational polynomial P{T) to each side of the equation: : 
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p^ryn^u ^ p(T)(^e"'^^y = (use 2.6) 

= e"'^*'P(e"''') = e"'^*'P((-l)") = P((-l)") cosriTrt + iP((-l)") sinnTri 
and 

Pijyn-^tt ^ p(T){cosmrt + ismmrt) = P{T) cosnirt + iP{T) smm^t 
Equating the real and imaginary parts and supposing that P(— 1) 7^ 0, we obtain 



1 COS H7rt 

(2.7) P{T) cosnirt = P {-!) cos nirt ^ cosiiTrt 



P(T) P((-l)") 



N 1 sinnTrt 

(2.8) P{T) smn-Kt = P{-1) smnnt — — sinnTrt 



P(T) P((-l)") 
Example 2.2. Find a particular solution yp{t) of the finite difference equation 
y{t + 2) - 5y{t + 1) + 6y{t) = cos(7ri) 

Solution: 

y{t + 2) - 5y{t + 1) + 6y{t) = cos(7rt) (T^ - 5T + 6)y(t) = cos(7rt) 

Replacing 2.7, with n = 1, we obtain 

, , 1 , , cos(7ri) , , cosint) 

= T^-5T + 6 = (-l)^-5.(-l) + 6 ^ = ^2- 

□ 

2.4. Equations where (/)(t) = A*f(t). Consider the polynomial 
P(T) = aoT" + aiT"-i + • • • + a„_iT + a„ 

We have 

P(T)A*/(i) = aoT"(A*/(t)) + aiT"-i(A*/(i)) + • • • + a„_ir(AV(i)) + a„AV(t) = 
= aoA*+"/(i + n) + aiA*+"- V(t + n - 1) + • • • + a„_iA*+V(t + 1) + anA*/W = 
= A*(aoA"T"/(t) + aiA"-ir"-V(<) + • • ■ + a„_iAr/(t) + a„/(t) = X'P{XT)f{t) 

Therefore, we obtain the property 



(2.9) P{T)X'f{t) = \'P{\T)f{t) ^ A*f (t) = A* (t) 
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Example 2.3. Find a particular solution yp{t) of the finite difference equation 
y{t + 2) - 5y(i + 1) + Ay{t) = 3* sin(7rt) . 

Solution: 

y{t + 2) - 5y{t + 1) + 4:y{t) = 3* sin(7rt) => (T^ - 5T + 4)y{t) = 3* sin(7rt) 
Replacing: 

1 

(3r)2 -T{3T) + 4 ' 



= ^ 9T^-15T + 4 ^^^('^^^ = ^ 9(-l)^-15(-l)+4 ^ ^^^^^ = ^Ss^ 

□ 

2.5. Formula for polynomials P(T — A) with ^(t) = A*f(t). Here I prove by 
induction that (T - A)"A*/(i) = A*[A(T - 

1. When n = 1 

iT-\)\\f{t) = \'+\f{t + l)-\'+\f{t) = A*+i(r-l)/(t) = A*[A(T-l)]/(t) 

2. Showing that the statement holds when n = k 

(r-A)'=A*/(i) = A*[A(r-l)]V(i) 

3. Prove when n = k + 1 

(T - A)'=+iA*/(i) = (T - A)(T - A)'=A*/(i) = (T - A)A*[A(r - = 

= X*+^[XiT-l)ff{t+l)-\'+^[\iT-l)ff{t) = A*+i[A(T-l)]'=[T/(t)-/(t)] = 

= A*+i[A(r - l)f{T - l)f{t) = A*+iA'=(T - if +^ fit) = A*[A(r - 
Finally, we generalicc for P{T — A)A*/(t). 

Consider the polynomial P{T) = aoT" + aiT""! H h a„_iT + a„. 

We have 

P(T - A)A*/W = 

= ao(r-A)"(A*/W)+ai(r-A)"-i(A*/W)+- • •+a„_i(T-A)(A*/W)+anA*/W = 
= aoA*[r(T-l)]"/(i)+aiA*[r(T-l)]"-V(i)+- ■ •+a„_iA*[r(T-l)]/(t)+a„A*/(i) 

= A*P[A(r-l)]/(t) 

Therefore, we can establish the property 

p(r - r)A*/W = A*P[A(r - ^ 

(2-10) <^ .x A*f(t) = A*— — -rf(t) 

^ ' P(T - A) ^ ^ P[A(T - 1)] ^ ^ 
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Example 2.4. Find a particular solution yp{t) of the finite difference equation 

y{t + 1) - 2y{t) = 2* 

Solution: 

y{t + 1) - 2y{t) = 2* ^ (T - 2)y{t) = 2* ^ y{t) = ^^2* 
We cannot apply 2.6, since it would give a division by zero. Instead, we write 

yW = T^(2*.l). 

By applying 2.10, we obtain 

= 2*^(1) = 2*1^(1) . 2-^(1) 

Finally, using 2.4: 

yp{t) = 2'-H 

□ 

3. Conclusions. This paper developed a technique for solving linear finite dif- 
ference equations with constant coefficients. In addition, it proves several fundamental 

properties (linearity, translation, unity and propagation) of the polynomial translation 
operator and establishes formulae for solving several forms of FDE: A* , cos nirt, sin nirt, A* 
and A*/(i). 
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